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ABSTRACT 


In  this  report,  we  discuss  adaptive  control  of  a  space  robot  system  with  an 
attitude  controlled  base  on  which  the  robot  is  attached.  We  at  first  derive  the 
system  kinematic  and  dynamic  equations  based  on  Lagrangian  dynamics  and  linear 
momentum  conservation  law.  Using  the  dynamic  model  developed,  we  discuss  the 
problem  of  linear  parameterization  in  terms  of  dynamic  parameters,  and  have  found 
that  in  joint  space  the  dynamics  can  be  linccirized  by  a  set  of  combined  dynamic 
parameters,  but  in  inertia  space  linear  parameterization  is  impossible  in  general. 
Then  we  propose  an  adaptive  control  scheme  in  joint  space  which  has  been  shown 
effective  and  feasible  for  the  cases  where  unknown  or  unmodeled  dynamics  must 
be  considered,  such  as  in  the  tasks  of  transport  an  unknown  payload,  or  catching 
a  moving  object.  The  scheme  avoids  the  use  of  joint  acceleration  measurement, 
inversion  of  inerti2d  matrix,  high  gain  feedback,  and  considerable  computation 
cost,  and  at  meamtime,  is  also  applicable  for  the  fixed-base  robot  system  by  slight 
modification.  Since  most  tasks  are  specified  in  inertia  space,  instead  of  joint  space, 
we  discuss  the  issues  associated  to  adaptive  control  in  inertia  space  and  identify  two 
potential  problems,  unavailability  of  joint  trajectory  since  mapping  from  inertia 
space  trajectory  is  dynamic  dependent  and  subject  to  uncertainty,  and  nonlinear 
parameterization  in  inertia  space.  We  approach  the  problem  by  making  use  of 
the  proposed  joint  space  adaptive  controller  and  updating  joint  trajectory  by  the 
estimated  dynamic  parameters  and  given  trajectory  in  inertia  space.  In  the  case 
study  of  a  planar  system,  the  linear  parameterization  problem  is  investigat  .  I,  the 
design  procedure  of  the  controller  is  illustrated,  and  the  validity  and  eff.::ctiveness 
of  the  proposed  control  scheme  is  demonstrated. 
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1  Introduction 


Considerable  research  efforts  have  been  directed  to  some  primary  functions  of  robots  in  space  appli¬ 
cations,  such  as  material  transport  [17],  simple  manipulation  [4],  basic  locomotion  [14],  inspection 
and  maintenamce  of  the  space  station  and  satellites  [1,  4].  The  adaptive  control  is  critical  for  the 
robot  system  subject  to  dynamic  uncertainty  in  these  tasks. 

For  material  transport  and  manipulation  tasks,  space  robots  have  to  face  uncertainty  on  the 
parameters  describing  the  dynamic  properties  of  the  grasped  load,  such  as  moments  of  inertia  or 
exact  position  of  the  mass  center.  In  most  cases,  these  parameters  are  unknown  and  thus  they 
can  not  be  specified  off-line,  in  inverse  dynamics  for  feedforward  compensation  or  any  model- 
based  control  scheme.  In  catching  a  moving  object  [15],  the  robot  is  expected  to  be  capable  of 
adaptation  to  the  dynamics  change  at  the  moment  of  catching  operation.  On  the  other  hand,  most 
space  robots  are  designed  to  be  light- weighted  and  thus  low-powered,  partially  due  to  the  fact 
zero  gravity  environment.  Therefore,  joint  friction  and  damping  in  space  robots  are  much  more 
significant  than  those  in  industrial  robots.  These  effects  are  neither  negligible  nor  easy  to  model. 
Adaptive  control  may  provide  a  feasible  solution  to  those  system  dynamics  uncertainties.  Adaptive 
control  is  also  able  to  accommodate  various  unmodeled  disturbances,  such  as  base  disturbance, 
microgravity  effect,  sensor  and  actuator  noise  due  to  extremes  of  temperature  and  glare,  or  impact 
effect  during  docking  or  rendezvous  process. 

Most  of  existing  adaptive  control  algorithms  have  the  following  shortcomings  which  cause  their 
applications  in  space  robots  unrealistic,  the  use  of  joint  acceleration  measurement,  the  need  of 
inversion  of  inertia  matrix,  high  gain  feedback  and  considerable  computational  cost.  The  first 
two  must  be  avoided  even  for  fixed-base  industrial  robots,  because  of  lack  of  joint  acceleration 
sensor  and  the  complexity  of  inversion  of  inertia  matrix.  Slotine  and  Li  [11]  have  tackled  these 
problems  successfully.  A  high  gain  feedback  is  extremely  harmful  for  a  space  robot  which  is  usually 
light- weighted  and  low-powered.  Considerable  computation  also  needs  to  be  avoided  for  allowable 
package  of  self-contained  space  robots. 

This  report  focuses  on  the  robot  system  where  the  base  attitude  is  controlled  by  either  thrust 
jets,  or  reaction  wheels.  The  reaction  wheels  are  arranged  in  orthogonal  directions,  and  the  number 
of  reaction  wheels  can  be  three  or  two  depending  on  different  tasks,  and  a  standard  reaction  wheel 
configuration  can  be  found  in  [6].  When  the  attitude  of  the  base  is  controlled,  the  orientation  and 
position  of  both  robot  and  base  are  no  longer  free,  and  the  dynamic  interaction  between  the  base 
and  robot  results  in  the  dynamic  dependent  kinematics,  i.e.,  the  kinematics  is  in  relation  to  the 
mass  property  of  the  base  and  robot.  Control  is  not  only  applied  to  robot  joint  angles,  but  also 
three  orientations  of  the  base. 

In  this  report,  based  on  linear  momentum  conservation  law  and  Lagrangian  dynamics,  we  at 
first  formulate  kinematics  and  dynamics  equations  of  the  space  robot  system  with  an  attitude 
controlled  base,  in  a  systematic  way.  Based  on  the  dynamic  model  developed,  we  study  the  linear 
parameterization  problem,  i.e.,  dynamics  can  be  linearly  expressed  in  terms  of  dynamic  parameters, 
such  as  mass  and  inertia.  We  have  found  that  for  the  space  robot  system  with  an  attitude  control 
base,  the  linear  parameterization  is  valid  in  joint  space,  while  is  not  valid  in  inertia  space  which 
can  be  viewed  as  Cartesian  space  for  earth-based  robots. 

Using  the  dynamic  model,  we  propose  an  adaptive  control  scheme  in  joint  space.  The  scheme 
does  not  need  to  measure  accelerations  in  joint  space,  and  a  high  feedback  gain  is  not  required.  The 
proposed  method  is  eflTective  and  feasible  for  space  robot  applications  when  dynamic  parameters 
are  unknown  or  unmodeled  dynamics  effect  must  be  considered.  Since  in  most  applications,  the 
tasks  are  specified  in  inertia  space  normally,  instead  of  joint  space,  we  discuss  the  issues  in  relation 
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Figure  1:  Space  robot  system  with  an  attitude  controlled  base. 

to  implementation  of  adaptive  control  in  inertia  space  and  identify  two  main  problems.  The  first 
problem  occurs  when  the  joint  adaptive  control  is  executed.  The  required  joint  trajectory  cannot 
be  generated  by  the  given  trajectory  in  inertia  space  due  to  the  parameter  uncertainty  in  the  kine¬ 
matic  mapping  which  is  dynamics  dependent.  The  second  problem  is  nonlinear  parameterization 
in  inertia  space  which  make  impossible  to  implement  the  same  structured  adaptive  control  as  that 
in  joint  space.  We  approach  this  problem  by  making  use  of  joint  space  adaptive  controller  and  up¬ 
dating  joint  trajectory  from  identified  kinematic  mapping  and  the  given  trajectory  in  inertia  space. 
This  method  has  shown  its  effectiveness  in  simulation,  and  some  issues  eissociated  to  parameter 
estimation  and  updating  time  are  discussed. 

Finally,  we  study  a  planar  robot  system  to  investigate  linear  parameterization  problem  of  robot 
system  dynamics,  smd  illustrate  the  validity  and  effectiveness  of  the  proposed  adaptive  control 
schemes. 

2  Kinematics  of  Space  Robot  System  with  An  Attitude  Con¬ 
trolled  Base 

In  this  section  we  discuss  the  kinematics  of  the  space  robot  system  when  the  orientation  of  the 
base  is  controlled  and  the  translation  of  the  base  is  free.  The  relationship  between  the  robot  hand 
motion  in  inertia  space  and  robot  joint  motion  is  derived  using  linear  momentum  conservation  law. 

As  shown  in  Figure  1,  a  space  robot  system  with  an  attitude  controlled  base  can  be  modeled 
as  a  multibody  chain  composed  of  n  -i-  1  rigid  bodies  connected  by  n  joints,  which  are  numbered 
from  1  to  n.  Each  body  is  numbered  from  0  to  n,  and  the  base  is  denoted  by  B  in  particular. 
The  mass  and  inertia  of  ith  body  are  denoted  by  mi  and  li,  respectively.  A  joint  variable  vector 
<1  =  (9ii92>'  ■  ■  is  used  to  represent  those  joint  displacements.  The  orientation  of  the  base  is 
represented  by  a  vector  qa  =  iqBi,QB2,qB3)^ ■ 

Two  coordinate  frames  are  defined,  the  inertia  coordinate  orbit,  and  the  base  co¬ 

ordinate  53b  attached  on  the  base  body  with  its  origin  at  the  centroid  of  the  base.  As  shown  in 
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Figure  1,  let  R,-  and  rj  be  the  position  vectors  pointing  the  centroid  of  ith  body  with  reference  to 
53/  and  5Ze  respectively,  then 


R,’  =  Pi  +  R-B  ( 1 ) 

where  Rb  is  the  position  vector  pointing  the  centroid  of  the  base  with  reference  to  53/-  I"®*- 
and  n,-  be  linear  and  angular  velocities  of  ith  body  with  respect  to  53/*  v,  and  u;,  with  respect  to 
53  b-  Then  we  have 


V,  =  V,  +  Vb  +  Qb  X  r,- 

fl,-  =  Ui  +  Qb  (2) 

where  Vb  and  Hb  are  linear  and  angular  velocities  of  the  centroid  of  the  base  with  respect  to  53/* 
and  operator  ’  x  ’  represents  outer  product  of  vector.  The  velocities  v,  and  w,  in  base  coordinates 
can  be  represented  by 


Vi  =  J£,i(q)q 


(3) 


Wi  =  J/4i(q)q 


(4) 


where  J£,i(q)  and  JMi(q)  are  the  submatrices  of  Jacobian  of  the  ith  body  representing  linear  part 
and  angular  part  respectively.  The  centroid  of  the  total  system  can  be  determined  by 


nc  =  ^rni 

i=0 

lc  =  th 
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(5) 

(6) 

(7) 

(8) 


The  linear  momentum  can  be  determined  by 

P  =  [Hv',Hvn] 


Vb 

fie 


+  Hv,q 


where 


=  HvVb  +  HvnfiB  +  Hv7q 

(9) 

Hv  =  nicUs 

(10) 

Hvn  =  -mc[rcx] 

(11) 

Hv,  =  mcJc 

(12) 
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and  U3  is  a  3  X  3  unity  matrix.  The  matrix  function  [rx]  for  a  vector  r  =  [ri,  ry,r*]^  is  defined  as 


[rx]  = 


0 


(13) 


Because  there  is  no  external  force  applied  to  the  system,  the  linear  momentum  is  conserved. 
However,  the  angular  momentum  is  not  conserved  for  attitude  control  torques  are  applied.  The 
linear  momentum  is  zero,  assuming  stationary  initial  condition. 


P  =  0  (14) 

Therefore,  we  may  represent  the  base  linear  velocities  by  base  angular  velocities  and  robot  joint 
velocities,  i.e.. 


Vs  =  -l/mc[Hvo,Hv,] 


SIb 

q 


(15) 


Now  we  derive  the  relationship  between  the  motion  rate  in  inertia  space  and  that  in  joint  space. 
For  position  control  tasks,  we  are  interested  in  controlling  three  orientations  of  the  base,  and  six 
generalized  displacements  of  the  robot  end-effector  simultaneously.  Control  actions  are  instead 
applied  at  n  robot  joints  and  three  base  attitudes.  We  therefore  define  V  and  6  as  generalized 
velocities  in  inertia  space  and  joint  space. 


V  =  [Os,  Vsp 


(16) 


^  =  ln8,q]^ 

where  Vs  is  the  velocity  of  the  robot  end-effector  in  inertia  space. 


(17) 


Vs  =  vs  +  Vs  -f-  Jls  X 

Since  the  velocity  of  the  end-effector  in  the  base  coordinates  is  determined  by 

Vfi  =  Jfiq 

where  Js  is  the  manipulator  Jacobian  with  respect  to  the  base  coordinates. 


Vs  =  J^q  -  1/mc  [Hvn.Hvg] 


Ob 

q 


-  [rsxjfls 


(18) 

(19) 


=  -  [r^x],  Js  - 


Ttlr 


Ob 

q 


~  [[(^c  ~  ~  ^c] 

Therefore,  the  motion  rate  relationship  between  joint  space  and  inertia  space  can  be 
introducing  a  special  Jacobian  matrix  N  which  differs  from  the  Jacobian  in  fixed-base 
generalized  Jacobian  in  a  completely  free- flying  space  robot  system. 


(20) 

obtained  by 
robot  or  the 


‘  ne 

Vs 


(21) 
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where 


and  O3  is  a  3  X  3  zero  matrix. 


N  = 


U3  O3 

J  rr  ^tE 


Jrr  —  [(l^c  *’e)x] 

JrE  =  —  Jc 


(22) 

(23) 

(24) 


3  Dynamics  of  Space  Robot  System  with  An  Attitude  Con¬ 
trolled  Base 

In  this  section  we  discuss  dynamics  of  the  space  robot  system  with  an  attitude  controlled  base. 
After  formulating  total  kinetic  energy  of  the  system  we  derive  the  dynamics  equation  of  the  system. 
Then  we  investigate  the  property  of  linear  parameterization  of  the  system  dynamics  which  is  critical 
for  developing  the  adaptive  control  algorithms  in  the  following  section. 

The  total  system  kinetic  energy  is  represented  by 


T  =  1/2Vb^HvVb  +  VB^lHvn.Hv,] 

+  l/2[n«,q|  f  1  f 


Hn 

Hn,  ■ 

Hb 

.  Hn,^ 

H, 

• 

q 

=  l/(2mc^)(nB,q]  Hv[Hvn,Hv,] 


-l/mo[nB,ql  [Hvn.Hv,] 

tlVg  q 


+  l/2[nB,q] 


Ho  Hn,  Hb 
Hn/  H,  q 


-  1 /ofo  •!  Hn  -  Hvn^Hvn/mc  Hn,  -  Hvn^Hv,/  1  [  Hb 
-l/2t«B,qJ  Hn,  -  Hvn^HK,/mc  H,  -  Hv/Hv,/m,  q 


=  l/2[nB,q]M(<>)  ^  J 

=  l/2e'^M{9)0  (25) 

where  M  is  the  inertia  matrix  of  the  system,  H,  is  the  robot  inertia  matrix  in  base  coordinate,  i.e., 
fixed  base  inertia  matrix,  and 

®  =  [qB,q]^  (26) 

n 

Hn  -  Ic  +  D(ri)m,-  (27) 
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(28) 


Hn,  =  +  »".[r.)JL.) 

t=i 


M($)  = 


Mil  Mi2 
M21  ^^22 


Mil  =  Ho  -  Hy^Hvo/mc 
Mi2  =  Hn^  —  Hy^Hyg/nic 
M21  =  Mji 

M22  =  H,  -  Hygilyg/mc 


(29) 

(30) 

(31) 

(32) 

(33) 


The  property  of  linear  parameterization  to  dynamic  parameters  is  one  of  prereqnisite  condi¬ 
tions  under  which  most  adaptive  and  nonlinear  dynamic  control  schemes  are  developed.  It  has  been 
shown  that  the  problem  of  parameterization  linearity  in  dynamics  can  be  reduced  to  the  problem 
of  parameterization  lineajity  in  inertia  matrix.  Therefore,  in  order  to  study  whether  linear  param¬ 
eterization  is  valid  for  the  space  robot  system  with  an  attitude  controlled  base,  we  need  to  show 
whether  the  inertia  matrix  M  can  be  linearly  represented  by  a  set  of  properly  chosen  combinations 
of  dynamic  parameters. 

Based  on  the  previous  derivation,  each  member  of  matrix  M  can  be  further  expanded  in  the 
following  forms. 
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^22  = 

i=l  f=l 

n  n 

■»33  =  (X!  +  (X 

i=l  i=l 

n  n 

Sl2  =  -(X  •  (X 

i=l  i=l 

n  n 

•»23  =  -(X  '"»’■»»)  •  (X 
1  =  1  1  =  1 

n  n 

«i3  =  -(X 


1=1 


«=i 


and 


(34) 

(35) 

(36) 

(37) 

(38) 

(39) 

(40) 
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D(r)  =  [rx]^[rx]  = 


rl  +  -TxTy  -r^r* 


-r-r 


l' V 


Thus, 


[  TyTt  Tj.  +  ry  J 

n  n  n 

Mil  =  Ic  +  H  D(r,)m.  - 

i=l  «=1  i=l 

n  n 

M22  =  H, - 


(41) 


(42) 


.=1 

n  n 


t=l 


.=1  j=i 


n  »*  **  ** 

Mi2  =  y\  +  53[rilJ£,.m.-  +  21^  SijTn,m,/m, 
i^i  .=1  t=i  j=i 


(43) 


(44) 


where  the  matrices  [r.]J/^.,  Ri,  Qi,  Si  are  only  functions  of  geometric  parameters,  i.e.,  indepen¬ 
dent  of  dynamic  parameters.  The  above  formulations  imply  that  the  inertia  matrix  can  be  linearly 
represented  by  a  set  of  combination  of  dynamic  parameters,  mjt,  Ifc,  mjmj/mc,  i,  j,  fc  =  0, 1,  •  •  • ,  n. 
From  the  kinetic  energy  formulation,  we  can  derive  dynamics  equation  by  Lagrangian  dynamics. 


M^4-B(^,0)^  =  r 


where 


^  1  oTi 


B{9j)e  =  M9-^i-0^Me) 

The  corresponding  dynamic  equation  in  inertia  space  is 

Hx  +  C(x,x)x  =  F 

where 

H  =  N"^MN"‘ 

C  =  N-^BN"‘  -HNN'* 


(45) 

(46) 

(47) 

(48) 

(49) 


N  is  a  generalized  Jacobian  matrix  and  is  dynamics  depenedent  for  the  space  robot  system.  The 
inertia  space  dynamic  equation  can  be  linearly  expressed  in  terms  of  dynamic  parameters  if  and 
only  if  the  inertia  matrix  H  can  be  linearly  parameterized  since 


C(x,x)x  =  Hx-|^(ix^Hx) 


(50) 


We  suppose  N"'  exists,  and 

jq-i  =  - 

det(N) 

where  N*  and  det(N)  are  the  adjoint  and  determinant  of  the  matrix  N,  then 


N*^MN* 

[det(N)]2 


(51) 


(52) 
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In  the  above  equation,  the  generalized  Jacobian  matrix  [det(N)]^  appears  as  the  denominator. 
From  derivation  procedure  of  the  N  in  the  last  section,  it  is  clear  that  the  N  is  time- varying  and 
highly  coupled  by  dynamic  parameters,  i.e.,  mass/inertia.  For  such  a  complicated  nonlinear,  time- 
varying  function  combining  with  dynamic  parameters  and  time- varying  joint  angles,  it  is  impossible 
that  every  element  of  N*^MN*  has  the  common  factor  [det(N)]*  at  every  instant. 

Even  if  the  above  statement  is  true,  there  is  still  possibility  to  linearly  parameterize  H,  provided 
that  the  numerator  can  be  linearly  parameterized  and  the  denominator  can  be  expressed  as  a 
product  of  two  scalar  functions  with  only  one  containing  dynamic  parameters,  i.e., 

d6f(N)  =  /i(m.-,/.)/2(/?,)  (.53) 

where  /i  is  a  function  of  dynamic  parameters  which  are  unknown  but  constant,  /2  is  a  function 
independent  of  any  dynamic  parameters.  This,  unfortunately,  is  impossible  in  general  due  to  high 
coupling  between  dynamic  parameters  and  joint  variables.  For  excimple,  two  DOF  generalized 
Jacobian  may  contain  the  following  simple  terms 

det(N)  =  mi3in(d\)  +  m2cos(02)  (54) 

Even  for  such  a  simple  form,  det(N)  cannot  be  decomposed  as  a  product  of  two  functions  with  one 
containing  mi  and  m2  only,  and  nor  can  [det(N)]*. 

This  may  raise  a  question  why  for  a  fixed-base  robot  the  similar  structured  adaptive  control 
can  be  implemented  in  Cartesian  space.  This  is  because  that  the  Jacobian  in  the  fixed-base  robots 
is  only  kinematic  dependent,  i.e.,  a  function  of  geometric  parameters  and  joint  angles.  Since  the 
dynamic  interaction  between  the  base  and  the  robot,  the  generalized  Jacobian  for  a  space  robot  with 
an  attitude  control  base  is  dynamics  dependent,  i.e.,  not  only  a  function  of  geometric  parameters 
and  joint  angles,  but  also  a  function  of  the  dynamic  parameters.  It  is  to  these  parameters  that 
we  aim  at  adapting  in  our  problem.  Therefore,  the  inertia  matrix  for  the  fixed-base  robot  can  be 
linearly  parameterized  for  dynamic  parameters  in  Cartesian  space,  while  for  a  space  robot  it  is 
impossible  in  inertia  space. 

Generally  speaking,  for  a  space  robot  with  attitude  controlled  base,  dynamics  can  be  linearly 
parameterized  in  terms  of  dynamic  parameters  in  joint  space,  but  it  cannot  in  inertia  space. 

4  Adaptive  Control  Scheme 

In  this  section,  based  on  the  dynamic  model  developed  and  the  property  of  parameterization 
linearity  in  Section  3,  we  discuss  the  adaptive  control  strategy  for  space  robot  system  with  an 
attitude  controlled  base. 

At  early  state,  adaptive  control  approaches  for  conventional  fixed-base  robot  manipulators 
are  based  on  unrealistic  assumptions  or  approximations  on  local  linearization,  time-invariant  and 
decoupled  dynamics  [3,  5].  These  assumptions  or  approximations  are  relaxed  after  some  results 
developed  in  the  context  of  parameter  estimation  (8).  Based  on  the  possibility  of  selecting  a 
proper  set  of  equivalent  parameters  such  that  the  manipulator  dynamics  depends  linearly  on  these 
parameters,  research  on  adaptive  robot  control  can  now  take  full  consideration  of  the  nonlinear, 
time-varying  and  coupled  robot  dynamics.  As  stated  in  [7],  all  three  kinds  of  adaptive  controllers 
in  use,  direct  [2,  11],  indirect  [10],  composite  adaptive  controllers  [12],  rely  on  the  possibility  of 
linear  parameterization  of  manipulator  dynamics. 

From  previous  discussion,  we  have  learned  that  the  dynamics  of  the  space  robot  system  in  joint 
space  is  linear  in  terms  of  a  set  of  combinations  of  dynamic  parameters.  Therefore,  this  set  of  new 
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combined  parameters  can  be  used  in  the  design  of  our  adaptive  controller.  This  leads  us  to  propose 
an  adaptive  control  algorithm  in  joint  space.  Since  a  unique  solution  may  be  found  from  inverse 
kinematics  of  the  robot  system  with  the  attitude  controlled  base,  adaptive  control  algorithm  in  joint 
space  is  feasible.  However,  this  is  not  true  for  a  complete  free-flying  space  robot  system. 

Let’s  recall  the  dynamic  equation  in  joint  space 


MB-bB{B,9)B=  T 

(55) 

We  define  a  composite  error  s 

s  =  Cp  -f-  CCp 

(56) 

Cp  =  Bj  -  B 

(57) 

ep  =  Bj-B 

(58) 

and  we 

also  define  modified  joint  velocity 

B'  =  B +  3 

(59) 

and  modified  joint  acceleration, 

ii  .0 

1 

(60) 

i.e.. 

B  =  B  +  8  +  3  =  {9d  —  Cp)  +  {(ip  Ccp)  -i-  (cp  -)-  Ccp) 

—  +  (C  +  l)ep  +  Cep  =  +  s  +  Cep 

(61) 

If  we  apply  the  following  control  law  in  joint  space, 

T  =  mb"  -1-  bb' 

(62) 

then 

MB  +  B{B,B)B  =  Mb"  +  BB' 

(63) 

i.e.. 

MB  =  -BiB,B)B  +  mb"  +  BB' 

(64) 

Defining  M  =  M  —  M,  B  =  B  —  B,  we  have 


Mcp  =  Mdj  -  Me 

=  M[e"  -  s  -  Cep]  -  [-B^  -H  mb"  -|-  B0'] 

=  M[d''  -  s  -  Cep]  -  [-B{9'  -  s)  +  mb"  B^'] 
=  -Mfl"  -  BB'  -  (M  -I-  B)s  -  MCe'p 
=  -YiB,lBj,Bj,Bj)a  ~  (M  +  B)s  -  MCop 

where 

Ya  =  MB"  -i-  BB' 

a  =  a  -  a 
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(65) 

(66) 


and  a  is  estimation  of  the  unknown  dynamic  parameters  of  the  space  robot  system  including  the 
robot,  the  base,  and  probably  the  payload  which  is  being  manipulated. 

We  now  design  our  adaptive  control  algorithm  using  Lyapunov  function  candidate 

V  =  l/2s^Ms  +  l/2a^ra 

where  the  matrix  F  is  diagonal  and  positive  definite.  This  yields 

V  =  1/28^M8  +  8^M8  4-  a^Fa 

=  1/28^M8  4-  8^M(ep  4-  C®p)  +  a^Fa 
=  -8^Ya  -  8^(M  4-  B)8  4-  1/2s^M8  4-  a^Fi 
=  -8^M8  4-  1/28^(M  -  2B)8  4-  a^(Fi  -  Y^s) 

If  we  use  adaptation  law 

i  =  F“*Y^s 

then 

V  =  -8^Ms  <  0  (69) 

« 

due  to  the  fact  that  the  matrix  M  -  2B  is  skew-symmetric,  and  M  is  positive  definite.  Therefore, 
the  system  is  stable  in  the  sense  of  Lyapunov,  because  Y  is  a  positive,  nonincreasing  function 
bounded  below  by  zero.  a(t)  and  a(t)  are  then  bounded,  and  s(t)  is  a  so-called  square  integrable 
or  £2  function  [13].  Provided  that  the  function  Y  is  bounded,  this  is  sufficient  for  the  purpose  of 
control  because  a(t)  converge  to  zero  as  the  £2  function  must  converge  to  zero  as  t  -*  00.  The 
parameter  estimation  error  a(t)  will  converge  to  zero  only  if  persistent  excited  input  is  utilized. 
The  output  error 

s  =  ep  +  Cep  (70) 

converges  to  zero,  which  in  turn  implies  that  Cp  — >  0  as  t  — >  00  since  C  is  positive.  We  can  now 
readily  conclude  our  adaptive  control  algorithm  in  Theorem  1. 

Theorem  1  For  the  dynamic  system  (55),  the  adaptive  control  law  defined  by  (62)  and  (68)  is 
globally  stable  and  guarantees  zero  steady  state  error  in  joint  space. 

The  composite  error  s  is  of  PD  type  structure  which  is  the  same  as  the  composite  error  defined  by 
Slotine  and  Li  [11].  In  general  the  PD  structure  control  adds  damping  to  the  system  but  the  steady- 
state  response  is  not  affected.  The  PI  structure  adds  damping  and  improve  the  steady-state  error  at 
the  same  time,  but  rising  time  and  settling  time  are  penalized.  To  improve  the  system  steady-state 
error,  in  the  proposed  2Mlaptive  control  algorithm,  the  PID  type  s  can  also  be  used.  Since  when 
the  PID  type  s  is  used,  the  order  and  type  of  the  system  is  increased  by  one,  the  steady-state  error 
is  decreased,  and  thus  the  system  is  more  robust  to  parameters  uncertainties  which  usually  cause 
a  significant  steady-state  error.  Moreover,  the  PID  type  s  allows  two  parameters,  instead  of  one. 
to  be  adjustable  in  order  to  achieve  a  desired  system  performance.  In  what  follows,  we  discuss  the 
stability  of  the  control  scheme  when  the  PID  type  s  is  employed. 

Define  ^ 

8  =  ep4-Ctep  +  C2  /  (71) 

Jto 
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Figure  2:  Block  diagram  of  adaptive  controller  in  joint  space. 


and  the  gains  C2  be  selected  such  that  the  eigenvalues  of  the  tracking  error  equation 

ep  +  Ci^p  +  CsCp  =  0  (72) 

have  negative  real  parts.  This  ensures  the  global  stability  of  the  system  when  s  converges  to  zero. 
Using  the  PID  type  s  and  the  same  definitions  of  6  and  6  ,  we  can  derive  that 

Mcp  =  Mdd  -  M9 

=  M[d"  -  8  -  CiCp  -  Czep]  -  (-B0  +  Mq  +  B^’] 

=  -Y{9,9,0Mk  -  (M  +  B)s  -  MCiCp  -  MCzCp 


where 

Ya  =  M9"  +  B9' 
a  =  a  —  a 

When  the  same  type  of  Lyapunov  function  is  used 

V  =  l/2s^Ms  +  l/2a^ra 


then, 


If  adaptation  law 
is  used,  then 


V  =  l/2s^Ms  +  s^Ms  +  a^Fa 
=  l/2s^Ms  +  s^M(ep  -h  CiSp  +  CzCp)  +  a^Fa 
=  -s^Ms  +  l/2s^(M  -  2B)8  +  a^(Ffi  -  Y^s) 

i  =  F-*Y^s 

V  =  -s^Ms  <  0 


(73) 

(74) 

(75) 


(76) 


(”) 


for  all  s  due  to  the  fact  that  the  matrix  M  —  2B  is  skew-symmetric,  and  C11  C2  >  0,  and  M  is 
positive  definite. 

A  block  diagram  of  the  proposed  control  algorithm  with  PD  type  s  is  shown  in  Figure  2.  Our 
adaptive  controller  is  conceptually  simple  and  easy  to  implement.  This  approach  does  not  require 
the  use  of  joint  accelerations  and  inversion  of  inertia  matrix.  Its  computational  cost  is  low  because 
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it  can  be  implemented  through  the  use  of  Newton-Euler  recursive  formulation.  It  can  be  seen  from 
Equation  (62)  which  has  the  same  structure  as  computed  torque  method,  that  the  control  law  can 
be  computed  efficiently  using  a  Newton-Euler  formulation  once  a  have  been  specified.  A  high  gain 
feedback  is  not  a  must  for  the  system  stability.  This  adaptive  approach  can  also  be  applied  to 
industrial  robot  control  by  a  slight  modification. 

5  Adaptive  Control  in  Inertia  Space 

In  this  section,  we  extend  our  joint  space  adaptive  control  approaches  to  the  problems  where  control 
variables  are  specified  in  inertia  space. 

Conceptually,  for  most  applications,  the  desired  robot  hand  trajectory  (i.e.,  position,  velocity 
and  acceleration)  must  be  specified  in  inertia  space.  For  example,  let’s  consider  catching  a  moving 
object  by  a  space  robot.  The  desired  trajectory  after  catching  depends  upon  the  tasks  and  the 
motion  trajectory  of  the  object  before  catching,  and  thus  must  be  specified  in  inertia  space.  In 
other  words,  as  in  the  case  of  fixed-base  robot  tasks  are  normally  specified  in  Cartesian  space,  tasks 
in  space  applications  are  unlikely  to  be  specified  in  joint  space.  Fortunately,  the  mapping  from 
robot  hand  position  in  inertia  space  to  displacements  in  joint  space  can  be  uniquely  determined  for 
space  robot  system  with  an  attitude  controlled  base,  which  differs  from  the  case  of  a  completely 
free-flying  space  robot  system.  This  unique  kinetic  relationship  has  been  first  studied  by  Longman 
et  al.  [9],  and  also  is  illustrated  by  a  planar  example  in  our  case  study. 

However,  the  unique  kinematics  relationship  can  only  be  determined  when  dynamic  parameters 
are  given,  for  this  relationship  is  indeed  dynamic  dependent.  When  some  dynamic  parameters  are 
unknown,  which  is  indeed  the  reason  why  we  come  to  adaptive  control,  the  mapping  cannot  be 
determined!  Therefore,  the  primary  difficulty  of  extending  our  approach  from  joint  space  to  inertia 
space  is  that  the  desired  trajectory  in  inertia  space  cannot  be  transformed  to  the  desired  trajectory 
in  joint  space  because  some  dynamic  parameters  are  unknown.  In  previous  discussion,  we  have 
utilized  a  desired  trajectory  in  joint  space,  as  other  researchers  have  done  [16],  without  giving  any 
explanation  about  how  the  trajectory  is  generated.  The  problem  is  not  significant  if  the  objective 
is  to  identify  dynamic  parameters,  but  is  important  if  the  objective  is  to  control  the  system. 

The  problem  can  be  resolved  if  the  same  structured  adaptive  control  scheme  can  be  implemented 
in  inertia  space.  This,  however,  is  not  feasible  because  the  proposed  adaptive  control  scheme  in 
joint  space  requires  that  the  dynamic  model  must  be  linearly  parameterized.  Therefore,  the  same 
type  of  the  control  scheme  cannot  be  developed  in  inertia  space.  As  has  been  known,  the  dynamic 
related  generalized  Jacobian  of  space  robot  makes  it  impossible  to  suitably  choose  a  set  of  dynamic 
parameters  such  that  the  inertia  space  system  dynamics  can  be  linearized.  That  is  why  the  same 
structured  adaptive  controller  in  joint  space  is  not  feasible  for  adaptive  control  in  inertia  space. 

We  approach  the  problem  in  the  following  way.  First,  given  trajectory  in  inertia  space,  we  use 
an  initial  estimation  of  dynamic  parameters  to  compute  initial  joint  trajectory.  Then  the  initial 
joint  trajectory  and  dynamic  parameters  are  used  in  the  proposed  joint  space  adaptive  control 
algorithms.  After  a  certain  period  of  time  we  update  the  system  dynamic  parameters  by  using 
new  estimated  ones  in  the  outer  loop  of  our  controller.  We  can  then  specify  more  precise  joint 
space  trajectory  based  on  these  new  parameters  and  the  inertia  space  trajectory.  Since  the  inertia 
space  trajectory  is  uniquely  determined  by  the  joint  space  trajectory  and  dynamic  parameters,  it 
can  be  shown  from  the  Jacobian  relationship  that  position  errors  in  inertia  space  converges  to  a 
given  boundary  if  position  errors  in  joint  space  and  parameter  errors  are  bounded,  provided  that 
the  robot  is  not  in  its  singularity  configuration.  The  control  scheme  is  illustrated  in  Figure  3. 

It  is  worthwhile  to  discuss  two  issues  in  the  implementation  of  the  proposed  control  scheme. 
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Figure  3:  Block  diagram  of  adaptive  control  scheme  in  inertia  space. 

First,  to  accurately  estimate  unknown  parameters,  a  persistent  excitation  (PE)  trajectory  is  re¬ 
quired  to  drive  the  robot  joints.  PE  trajectories  in  joint  space  and  in  inertia  space  are  not  equivalent 
because  the  spectrum  of  trajectory  signal  in  inertia  space  is  different  from  the  spectrum  of  the  same 
signal  in  joint  space  due  to  nonlinear  kinematic  transformation.  Therefore,  it  is  of  importance  to 
carefully  choose  initial  trajectory  in  inertia  space  such  that  the  same  trajectory  in  joint  space  is 
PE.  If  the  PE  condition  is  not  satisfied,  parameter  identification  error  occurs,  although  the  joint 
space  position  errors  may  stiff  converge. 

Second,  the  updating  time  for  inverse  kinematics  using  the  estimated  parameters  in  outer  loop 
of  our  controller  must  be  slow  enough  to  maintain  the  system  stable.  The  outer  loop,  as  shown  in 
Figure  3,  is  used  to  update  the  inverse  kinematics  and  therefore  the  desired  joint  trajectory  which 
is  used  in  joint  space  adaptive  controller.  A  fast  updation,  especially  using  incorrect  parameters 
Pi,  may  not  guarantee  the  convergence  of  parameter  errors.  In  the  simulation,  the  updating  time 
for  inverse  kinematics  is  set  to  10  seconds.  Simulation  results  have  shown  that  position  errors  in 
inertia  space  converge  to  zero  as  errors  in  joint  space  converge  to  zero  and  estimated  parameters 
converge  to  their  true  values. 

In  fact,  if  the  updating  time  for  inverse  kinematics  is  long  enough,  we  can  also  view  the  control 
scheme  as  a  two-phase  approach,  parameter  identification  phase  and  control  phase.  That  is,  to 
estimate  dynamic  parameters  in  joint  space  using  the  joint  space  trajectory  transformed  by  the 
given  inertia  space  trajectory  and  initial  guess  of  parameters,  then  to  control  the  system  in  inertia 
space,  once  the  dynamic  parameters  has  been  correctly  identified.  If  the  dynamic  parameters  are 
estimated  ideally,  the  control  phase  may  also  be  executed  using  dynamic  control  algorithm. 


6  Simulation  Study 

In  previous  discussion,  we  studied  kinematics  and  dynamics,  and  presented  an  adaptive  algorithm 
in  joint  space  for  a  general  multiple-degrees-of-freedom  space  robot  system  with  an  attitude  con¬ 
trolled  base.  In  this  section,  we  conduct  a  case  study  to  show  the  computation  of  the  proposed 
algorithms  and  their  feasibility  in  robot  motion  control.  Though  the  following  discussion  is  con¬ 
fined  to  adaptation  to  mass  variation  only,  our  algorithm  is  also  applicable  to  other  parameter 
adaptation,  provided  that  a  set  of  combinations  of  those  parameters  can  be  chosed  such  that  the 
dynamics  can  be  linearly  expressed  in  terms  of  the  parameters  interested. 

A  two-DOF  revolute  manipulator  with  link  length  given  by  li  and  /j  {li=l2=l)  is  considered 
as  a  lumped-parameter  model  with  point  mass  mi  and  mj  at  the  end  of  each  link.  For  simplicity, 
we  assume  that  the  base  attitude  can  be  successfully  controlled  so  that  we  need  only  consider  the 


Figure  4:  A  planar  space  robot  system  model. 


control  of  the  robot  itself.  However,  it  must  be  pointed  out  that  our  adaptive  control  algorithm 
can  be  applied  to  control  the  robot  motion  and  the  base  orientation  simultaneously,  albeit  more 
complicated.  The  system  model  for  simulation  study  is  shown  in  Figure  4. 

At  initialization,  me  and  Re  are  computed,  and  they  remain  unchanged  unless  a  load  is  added. 


me  =  mo  +  mi  +  m2 

(78) 

mcRc  =  moRo  +  miR|  +  m2R2 

(79) 

Ri  =  Ro  +  ri 

(80) 

R2  =  Ro  +  r2 

(81) 

When  the  robot  is  in  motion. 

Ro  =  Rc  -  — ri  -  — r2 

(82) 

TTXq 

Re  =  R2 

(83) 

The  generalized  Jacobian  is 

N  =  Jfi  -Je 

(84) 

and 

Je  =  J2 

(85) 

(86) 

TTlc 

l  ~(*Tii  +  Tn2)s\  —  m2Si2  — m2Si2 

me  (mi  4-  m2)ci  +  m2Ci2  m2Ci2 

(87) 

where  s  and  c  stand  for  sine  and  cosine,  e.g.,  Sj  =  sin(qi),  ci2  =  cos{qi  +92)-  The  system  dynamics 
has  the  following  form, 

Mq  +  B(q,q)q=T  (88) 
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where 


therefore, 


TTlr 


/2 

M  =  — 

TUr 


M  =  M,  -  Ma 

„  _  .2  nii  +  2m2(l  +  C2)  m2(l+C2) 

’  "12(1 +  C2)  m2 

M2  =  mcJ^Jc 

TTij  -f-  2m2(mi  +  m2)(l  +  ca)  m2(l  +  ca)  +  mimaca 


m2(l  +  Ca)  +  TrtimaCa 


mi, 


momi  +  mima  +  2mom2(l  +  ca)  mi  m2  +  mom2(l  +  ca) 
mima  +  momaCl  +  Ca)  mima  +  mamo 

=  PiRi  +  P2R2  +  PsR-a 


where 


Ri  = 


0 

0  0 


Ra 


Pi  = 

nic 

P2  = 

mi  m2 

'Tic 

P3  = 

mo  m2 

me 

■ 

P  ' 

R3  = 

/2 

P 

2(1  +  C2)/2  (1+C2)/2 

(1+C2)/2 


(89) 

(90) 

(91) 

(92) 

(93) 

(94) 

(95) 

(96) 

(97) 

(98) 


It  is  noted  that  M  is  linear  in  terms  of  combined  dynamic  parameters  pi,p2  and  p3.  This  is  an 
example  to  show  that  dynamics  of  the  space  robot  system  with  an  attitude  controlled  base  can  be 
linearly  parameterized  in  joint  space.  We  also  note  that  mo,  mi  and  m2  can  be  uniquely  determined 
by  pi,p2  and  pa, 

.1  1  1  , 

mi  =  PiM—  +  —  +  — ) 

P\  P2  P3 


m2  =  P2P3(—  +  —  +  — ) 
Pi  P2  P3 

,1  1  1  , 

mo  =  PiP3(—  +  —  +  — ) 
Pi  P2  P3 


The  matrix  B  is  determined  by 

mo  m2 


B  = 


rric 


—  2/^3292  -^*^292 

l^S2qi  0 


=  P3R4 


where 


—  2/*S2?2  ~^^^292 

/*32?i  0 


Our  adaptive  control  law  is 


R4  = 

r  =  Mq  +  Bq’  =  Ya 


Y  = 


Riq  R2q  Raq  +  R4<1 
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(99) 

(100) 

(101) 

(102) 

(103) 

(104) 

(105) 


with  the  following  adaptation  law 


Pi 

a  =  P2 
P3 


7iS^R,q" 
728^R2q" 
738^(R3q"  +  R4q') 


(106) 


(107) 


To  study  the  proposed  adaptive  algorithms,  we  use  the  following  common  set  of  conditions: 

qu  =  +  6(8»n(0  +  co8(4t)))  ( 108) 

92d  =  :^(-126  +  6(stn(20  +  cos(60))  (109) 

C  =  10  (110) 

In  the  first  case  we  used  the  following  mass  parameters,  mo  =  41fcj,  mi  =  5kg,  m2  =  ‘ikg, 
and  the  initial  guess  of  aU  three  parameters  is  set  to  50%  of  their  true  values.  It  can  be  found 
from  Figure  5  that  joint  errors  converge  to  zero  and  all  parameters  converge  to  their  true  values 
4.1,  0.4,  and  3.28  (with  small  relative  errors  1.2%,  2.1%,  2.5%,  respectively)  after  a  transient 
period  (approximately  10  seconds).  The  results  showed  the  validity  and  efficiency  of  the  adaptive 
algorithm  proposed. 

We  then  compare  the  performance  of  adaptive  controller  and  dynamic  controller  without  adap¬ 
tation  when  there  is  uncertainty  in  dynamic  parameters.  In  order  to  make  the  dynamic  control 
more  favorable,  we  use  80%  of  true  values  as  initial  estimates  of  those  dynamic  parameters.  The 
dynamic  control  algorithm  is  based  on  PD  type  structure  in  joint  space  without  consideration  of 
parameter  uncertainty.  Figure  6  gives  plots  of  the  variations  of  two  joint  position  errors  by  using 
adaptive  control  and  dynamic  control.  The  adaptive  control  performance  is  distinctly  superior  to 
the  dynamic  control  response. 

To  study  the  effect  of  mass  ratio  of  the  base  with  respect  to  the  robot,  we  performed  simulation 
when  the  base  mass  is  sufficiently  large  compared  to  that  of  robot.  Figure  7  gives  the  simulation 
results  when  the  base  mass  is  50000kg.  The  results  have  shown  that  the  performance  is  not  sensitive 
to  the  mass  ratio,  which  also  shows  that  the  proposed  control  algorithm  is  applicable  to  fixed- base 
robots. 

Figure  8  shows  identification  of  combined  parameters  p\,  p2,  and  p3,  and  the  resultant  mass 
mj,  m2,  and  mo,  in  the  above  case.  From  Figure  8  we  found  that  estimation  of  all  parameters  mi, 
m2,  and  mo  are  very  close  to  their  true  values.  This  demonstrated  that  identification  of  combined 
dynamic  parameters  is  equivalent  to  the  identification  of  dynamic  parameters  mi,  m2,  and  mo,  as 
we  have  discussed  previously.  It  is  interesting  to  note  that  in  Figure  8  the  estimation  of  nonlinear 
dynamic  parameters  pi,  and  p^  converged  to  mg  and  m2  due  to  the  fact  that  the  base  mass  is 
almost  infinite. 

In  order  to  compare  two  different  adaptive  control  algorithms,  PD  type  and  PID  type,  various 
cases  have  been  tested.  For  a  persistent  excitation  (PE)  trajectory,  both  algorithms  presented 
almost  identical  performance.  For  a  non-PE  trajectory,  such  as 


qu  =  -  t  +  0.05t^) 
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92<i  =  j^(-120-^  +  0.05<2)  (112) 

the  steady  state  performance  is  improved  signihcamtly  using  PID  type  adaptive  controller,  as  shown 
in  Figure  9. 

For  inertia  space  adaptive  controller,  an  initial  guess  of  the  updating  parameters  is  set  to 
80%  of  the  true  value.  The  inertia  space  trajectory  and  joint  space  trajectory  employed  in  the 
simulation  are  shown  in  Figure  10.  We  used  10  seconds  as  updating  time  for  inverse  kinematics. 
The  effectiveness  of  this  adaptive  scheme  has  been  verified  by  the  tracking  errors  shown  in  Figure 
11.  It  is  found  that  position  errors  in  inertia  space  converge  to  zero  as  errors  in  joint  space  converge 
to  zero  and  estimated  parameters  converge  to  their  true  values. 

7  Conclusions 

In  this  report,  we  have  discussed  adaptive  control  of  a  space  robot  system  with  an  attitude  controlled 
base  on  which  the  robot  is  attached.  Adaptive  control  is  critical  for  various  robotic  applications  in 
space,  such  as  material  transport  and  light  manipulation,  in  which  robots  have  to  face  uncertainty 
on  the  dynamic  parameters  of  the  load  or  the  structure.  Based  on  Lagrangian  dynamics  and  linear 
momentum  conservation  law,  we  derived  system  dynamic  equations.  Then  we  showed  that  the 
system  dynamics  in  joint  space  can  be  linearly  parameterized,  i.e.,  the  dynamics  can  be  linearized 
in  joint  space  by  a  set  of  combined  dynamic  parameters,  while  the  same  conclusion  is  not  true  in 
inertia  space. 

An  adaptive  control  scheme  in  joint  space  is  proposed  to  cope  with  dynamic  uncertainties  based 
on  the  dynamic  model  developed.  The  scheme  is  effective  and  feasible  for  space  robot  applications 
by  eliminating  the  use  of  joint  acceleration  measurement,  inversion  of  inertial  matrix,  high  gain 
feedback,  and  considerable  computation  cost.  At  meantime,  the  scheme  is  also  applicable  for  the 
fixed-base  robot  system  by  slight  modification. 

Considering  that  the  tasks  in  space  are  specified  in  inertia  space  in  most  applications,  we 
discussed  the  issues  of  adaptive  control  of  the  robot  for  the  tasks  that  must  be  filfuU  in  inertia 
space.  Two  main  problems  have  been  identified.  If  the  joint  adaptive  control  is  implemented,  the 
desired  joint  trajectory  cannot  be  generated  from  the  given  inertia  space  trajectory  since  kinematic 
mapping  is  dynamics  dependent,  and  thus  is  subjected  to  uncertainty  in  parameters.  Moreover,  the 
same  structured  adaptive  control  as  in  joint  space  is  not  feasible  for  inertia  space  due  to  nonlinear 
parameterization  in  inertia  space.  We  approached  this  problem  by  making  use  of  the  proposed 
joint  space  adaptive  controller  while  updating  joint  trajectory  by  using  the  estimated  dynamic 
parameters  and  the  given  trajectory  in  inertia  spaK;e.  This  method  has  shown  its  effectiveness  in 
simulation.  Parameter  estimation  and  updating  time  are  discussed. 

Finally,  a  planar  system  is  studied  numerically  to  investigate  the  linear  parameterization  prob¬ 
lem  and  illustrate  the  procedure  to  design  the  controller.  The  results  demonstrated  validity  and 
effectiveness  of  the  proposed  adaptive  control  schemes  in  both  joint  and  inertia  space  descriptions. 
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Figure  5 
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Figure  7  Example  of  adaptive  control  for  fixed- based  robot 
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Figure  8  Illustration  of  combined  dynami 
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Trajectories  in  joint  space  and  in  inertia  space 
using  inertia  space  adaptive  controller 


Figure  11  Tracking  errors  in  joint  space  and  in  inertia  space 
using  inertia  space  adaptive  controller 
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